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Abstract. For an ideal D-dimensional Fermi gas under generic external confinement 
we derive the correcting coefficient (D — 2)/3D of the von Weizsacker term in the 
kinetic energy density. To obtain this coefficient we use the Kirzhnits scmiclassical 
expansion of the number operator up to the second order in the Planck constant h. 
Within this simple and direct approach we determine the differential equation of the 
density profile and the density functional of the Fermi gas. In the case D — 2 we find 
that the Kirzhnits gradient corrections vanish to all order in h. 
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1. Introduction 



At ultra-low temperatures macroscopic effects of quantum statistics are observable with 
bosonic vapors of neutral alkali-metal atoms [H[2]. In the last years quantum degeneracy 
has been achieved also with fermionic atoms [31 SI [5l El [7j. For fermions, it has been 
predicted that a reduced dimensionality strongly modifies density profiles [B~l l9l flOl ITT], 
collective modes [12] and stability of mixtures [131 El- I n a U these papers a Thomas- 
Fermi- von Weizsacker density functional approach is used [T5l fT6]. 

Some years ago, Holas, Kozlowski and March [17] derived for a D-dimensional ideal 
Fermi gas a correcting coefficient to the von Weizsacker term in the kinetic energy 
density. They used a response-function approach based on the knowledge of the static 
susceptibility in a D-dimensional space [18]. The aim of this paper is to recover the 
result of p2] by using the Kirzhnits semiclassical expansion of the number operator 
[T9] . We shall show that the Kirzhnits expansion in powers of h is a direct and simple 
approach to get the density profile and the energy of a D-dimensional ideal Fermi gas 
under generic external confinement. 



2. Thomas-Fermi- von Weizsacker formulation 



For an ideal D-dimensional spin-polarized Fermi gas in a hypercubic box of side L the 
number of particles at zero temperature is given by 
Ar rd D rd D p^{ p 2 \ 

N= - — -f 9 u-— , 1) 

where O(x) is the Heaviside step function. The chemical potential fl, that is the Fermi 
energy, fixes the number of particles N. In the formula we have substituted the discrete 
summation over box quantum numbers with the 2D-dimensional integral over single- 
particle phase-space variables (r,p). This approximation becomes exact in the limit 
L, N — > oo and n = N/L D constant. 

One can easily integrate over r and p in Eq. (CQ) and finds the following expession 
for the uniform density 

D/2 D 

2T(f + l) 
where 

/■oo O 

Fd{P) = J q dyy D ^Q^-y) = ^ D l\ (3) 

measures the number of states with energy less or equal to \x. 

In many papers the inclusion of an external potential U(r) is treated within the 
so-called local density approximation, i.e. by simply using a local chemical potential 

Mr) = /i - U(r) (4) 

in the Eq. ([5]), which becomes 



n 



/ m x D/2 D ^ . . 
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This expression of the local density n(r) can be obviously obtained from Eq. ([I]) by 
inserting on it Eq. and integrating over momenta p. Eq. (jSJ), also called Thomas- 
Fermi formula of the density profile [15], can be written in the following form 

9 *2 , /-) N 2/D 

+ l) n(r) 2 ^ + C/(r) = fi , (6) 
m V 2 J 

to make explicit the dependence of the chemical potential ft on the local density n(r). 
It is important to stress that Eq. (j5J), and also Eq. (jSJ), is semiclassical because 
it is derived by neglecting the fact that in quantum mechanics the kinetic operator 
T = —h 2 V 2 /(2m) and the operator fi = /i(r) — ft, — U(r) do not commute. Historically, 
it was von Weizsacker [16] the first to introduce a gradient correction (in the case D = 3) 
which takes into account phenomenologically the increase of kinetic energy due spatial 
variation of the density. Including the gradient correction of von Wiezsacker in Eq. ([6]) 
one has 



ri(Vn(r)) 2 lV 2 n(r) 

2m 



f r - + 1 n{vf' D +U{v) = fi. (7) 

m V 2 J 



4 n(r) 2 2 n(r) 

Remarkably, the gradient correction of von Weizsacker does not depend on the 
dimensionality of the system, as shown in [T7] . 



3. Kirzhnits gradient expansion 

Let us now consider a full quantum- mechanical treatment. In quantum mechanics the 
number operator is 

N = e(ji-f) . (8) 

The exact local density n(r) can be written by using the number operator N, the 
eigenstate |r) of the position operator, and the eigenstate |p) of the momentum operator 
in the following way 

n(r) = <r|iV» = J 7^^^^'^ . ( 9 ) 
where the last expression has been obtained by inserting the completeness 

/A>Ip>(pI = i, (io) 

and the formula (p|r) = e~ ip ' r / h /(2%h) D / 2 . According to Kirzhnits [TU [20], one can 
expand the operator N at the second order of h in the following way 

n\p) = e (a - f ) | P > = e(A - t) | P > + \&"{fL - t) [a, f ] | P ) 
- - t) [[ft, [ft,f]) + [f , [fc,f]]} | P > + l -&"\fi - m,f] 2 \ P ) (ii) 

where [fi, T] = fYT—Tfi and T\p) =p 2 /{2m)\p) =T\p). Here Q'(x) is the first derivative 
of Q(x) with respect to x, i.e. Q'(x) = 5(x) is the Dirac delta function, Q"(x) is the 
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second derivative of Q(x), and so on. It is important to observe that Eq. (11 ip does not 
depend on the dimensionality of the system. Kirzhnits used Eq. (fl~Tj) in the special case 
D = 3 [19J, while here we consider it for a generic D. Note that higher order terms 
of the Kirzhnits expansion of the number operator N involve higher derivatives of the 
Heaviside function 6(/t — T) PI ED]. 

The term (r|iV|p) in Eq. (Q can be now calculated by using the Kirzhnits 
expansion. In this way, after straightforward but lengthy manipulations, at the second 
order in h Eq. ([9]) becomes 

" (r) = (aS?) 2f(fTT){ FD( ' ,(r)) + 2^6 v VW^SWr)) 

+^(VMr)) 2 F£(Mr))]}, (12) 

where F' D {fi) is the first derivative of Fd(/j,) with respect to //, F^(fi) is the second 
derivative, and so on. In the Kirzhnits expansion, after integration over momenta, each 
term containing a n-th derivative of the Heaviside function 0(/2 — T) gives rise to a 
term containing a (n — l)-th derivative of the function Fd(/a). As a consequence, in the 
special case D = 2, where F 2 (fi) = fi, only the leading term of the Kirzhnits expansion 
is different from zero, i.e. with D = 2 the Kirzhnits gradient corrections vanish to all 
order in %. 

In the case D = 3, i.e. the three-dimensional Fermi gas, Eq. ( fl~2l) becomes 



_ p^(r)\ 3/2 J_ ( VMr) _ 1 (VMr)) 2 i . f . 

that is exactly the formula reported in the book of Dreizler and Gross [20J, which 
considers only the three-dimensional case. In the case D = 2, i.e. a two-dimensional 
Fermi gas, we have seen that all the quantum corrections are zero, while in the D = 1 
case these corrections are in general finite, and in particular the h 2 ones. 

We observe that the vanishing of quantum corrections in the Kirzhnits expansion for 
the two-dimensional ideal Fermi gas is a consequence of the linear number of states, i.e. 
constant density of states, in this specific dimension. Obviously, this does not mean that 
with D = 2 the Thomas- Fermi formula of the density profile is exact; it simply means 
that the Kirzhnits expansion fails in the D = 2 case. For instance, recently Brack 
and van Zyl [21] have found that, for the a 2D ideal gas of fermions under hamonic 
confinement, the quantum kinetic energy is exactly reproduced by the Thomas- Fermi 
kinetic energy without gradient corrections, but only using the exact density of states, 
which is different from the Thomas- Fermi one |21j . 

Eq. (fT2l) has a serious turning point problem: The local chemical potential 
/i(r) = ft—U(r) becomes zero at the classical turning points, where ft = U(r). Therefore, 
the gradient terms in Eq. ffl2|) diverge at the turning points. To overcome this problem 
we invert Eq. ( [121) finding /x(r) as a function of n(r). Obviously, due to the gradients 
this inversion cannot be done exactly but only in a perturbative way. First we note that 
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Eq. (1T21) has the following structure 

n(r) = A/,(r)^ 2 (l + I^BMr), VAi(r)]) , (14) 
where A = (m/(27rft 2 )) D / 2 /r(D/2 + 1). Then from this expression we get 

A -2/D n(r) 2/D 



Mr) 



^,V/i(r)] 

r 2 h 2 



2/D 



= A- 2 /^(r) 2 / D |l--— 5[Mr),VMr)] + 0(^ 4 )| , (15) 

where 0(h 4 ) contains terms of order h 4 . Finally, we find, at the second order in h, this 
expression 

h 2 (D-2\ Tl(Vn(r)) 2 lV 2 n(r)' 
2m" V 3D ) 4 n(r) 2 2 n(r) 

9 *2 , n N 2/D 

-" ft . r (- + l) n (r) 2 / D + U(r)= fx. (16) 



m V 2 

This result differs from that of von Weizsacker, Eq. (jTJ), simply because there is a the 
factor (D — 2) /3-D in the gradient terms. We stress that the von Weizsacker gradient 
terms are purely phenomenological, while the Kirzhnits gradient terms are based on a 
perturbative expansion of the exact number operator. In the scientific literature often 
the von Weizsacker gradient terms are modified by inserting a coefficient c in front of 
them. For instance, with D = 3, March and Tosi [22] suggest c = 1/3, while Zaremba 
and Tso [23] suggest c = 1/36. Remarkably, both the phenomenological improvements 
disagree with the Kirzhnits result c = 1/9. 

As previously discussed, the coefficient (D — 2) /ZD was found some years ago by 
Holas, Kozlowski and March [17] by using a response function approach based on the 
knowledge of the static susceptibility in a D-dimensional space |18) . Here we have 
instead applied the more direct Kirzhnits expansion. For D = 3 the factor gives the 
familiar coefficient 1/9, for D = 2 the factor is zero, and for D = 1 the factor is —1/3. 

It is straightforward to show that the differential equation (ITBl is obtained by 
minimizing the energy density functional 

]D r (D-2\ h 2 (Vn(r)) 2 



E[n(r)} = J d D r{[- 



+ 



3D / 8m n(r) 
27r^ 2 D /D \ 2 /° 



r(- + l) n(r) (D+2)/D + t/(r)n(r)} (17) 



m D 
with the constraint 

N = I d D rn(r) , (1* 



which fixes the chemical potential /i. Eq. (jTTI) can be used as a starting point to 
investigate mixtures of interacting Fermi atoms and also electrons in nanotubes. 
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4. Conclusions 

We have been able to determine the equations for the energy and the density profile 
of a D- dimensional ideal Fermi gas with external confinement by using the Kirzhnits 
expansion. These equations, which are correct at the second order in the Planck constant 
H, can be improved by including higher order terms of the Kirzhnits expansion of the 
number operator [19]. In fact, in some cases, ^-expansions can be calculated to all orders 
and resummed to give the exact result [21]. In the special case D = 2 we have found, 
by using the Kirzhnits expansion, that all h corrections to the Thomas-Fermi formula 
of the density profile are zero. This is a direct consequence of the constant density of 
states in a two-dimensional ideal gas and it means that the Kirzhnits expansion is not 
reliable under these conditions. We remind that gradient corrections to the Thomas- 
Fermi formula of the density profile have been analyzed with different semiclassical 
techniques in previous studies [T71 [25] . On the other hand, to our knowledge, this is the 
first time that the leading terms of the Kirzhnits expansion have been directly obtained 
for an ideal Fermi gas in a D-dimensional space and generic confining potential. 
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